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Abstract 

Background: Glutamyl queuosine-tRNA Asp synthetase (GluQ-RS) is a paralog of the catalytic domain of 
glutamyl-tRNA synthetase and catalyzes the formation of glutamyl-queuosine on the wobble position of tRNA Asp 
Here we analyze the transcription of its gene in Shigella flexneri, where it is found downstream of dksA, which 
encodes a transcriptional regulator involved in stress responses. 

Results: The genomic organization, dksA-gluQ-rs, is conserved in more than 40 bacterial species. RT-PCR assays 
show co-transcription of both genes without a significant change in transcript levels during growth of S. flexneri. 
However, mRNA levels of the intergenic region changed during growth, increasing at stationary phase, indicating 
an additional level of control over the expression of gluQ-rs gene. Transcriptional fusions with lacZ as a reporter 
gene only produced (3-galactosidase activity when the constructs included the dksA promoter, indicating that 
gluQ-rs do not have a separate promoter. Using bioinformatics, we identified a putative transcriptional terminator 
between dksA and gluQ-rs. Deletion or alteration of the predicted terminator resulted in increased expression of 
the lacZ reporter compared with cells containing the wild type terminator sequence. Analysis of the phenotype 
of a gluQ-rs mutant suggested that it may play a role in some stress responses, since growth of the mutant was 
impaired in the presence of osmolytes. 

Conclusions: The results presented here, show that the expression of gluQ-rs depends on the dksA promoter, 
and strongly suggest the presence and the functionality of a transcriptional terminator regulating its expression. 
Also, the results indicate a link between glutamyl-queuosine synthesis and stress response in Shigella flexneri. 

Keywords: tRNA modification, Gene expression, Stringent response, Osmotic stress 



Background 

The fidelity of the translation process depends on the 
aminoacyl-tRNA synthetase enzymes (aaRS). These es- 
sential enzymes are responsible for the correct attach- 
ment of the corresponding amino acid onto the cognate 
tRNA, therefore organisms have at least 20 synthetases 
[1]. The enzymes are divided in two classes, each class 
having a conserved structure. The genes encoding the 
aaRS are easily detected within sequenced genomes 
[2,3], and some species contain synthetase gene 
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duplications, such as the glutamyl-tRNA synthetases 
(GluRS) in Acidithiobacillus ferrooxidans and Helicobac- 
ter pylori (genes gltXl and gltX2) [4,5]. aaRS paralogs, 
predicted sequences with homology to fragments of 
synthetases, have also been identified, which is not sur- 
prising given the modular nature of the aaRS [6]. Some 
of the paralogs may be pseudogenes while others have 
known functions. For instance HisZ from Lactococcus 
lactis, which has similarity with the catalytic domain of 
histidyl-tRNA synthetase, is involved in histidine bio- 
synthesis [7]. A recent study in Salmonella enterica 
has shown that PoxA, encoded by poxAlgenX, has simi- 
larity to the carboxy-terminal catalytic domain of lysine- 
tRNA synthetase and is required for posttranslational 
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aminoacylation of bacterial elongation factor P. A poxA 
mutant has reduced colonization and virulence, possibly 
due to misregulated expression of proteins encoded by 
the SPI-1 pathogenicity island [8,9]. 

An Escherichia coli glutamyl-tRNA synthetase paralog, 
glutamyl queuosine-tRNA Asp synthetase (GluQ-RS) has 
approximately 35% amino acid similarity with the cata- 
lytic domain of GluRS. This includes the amino acids 
involved in recognition and activation of glutamate. 
Although GluQ-RS is missing the carboxyl-terminus do- 
main responsible for the tRNA recognition, in E. coli this 
enzyme is able to activate the amino acid in the absence 
of the tRNA. Further, once the aminoacyl- adenylate has 
been formed, the enzyme attaches the glutamate to the 
nucleoside queuosine present onto the tRNA Asp . There- 
fore, this enzyme is involved in the synthesis of a new 
modified nucleoside glutamyl-queuosine (GluQ) present 
in tRNA Asp [10,11]. This modification is present in 
tRNA isolated under acidic conditions from bacterial 
cells grown in rich media. However, the enzyme is not 
essential for growth of E. coli in rich or minimal media 
[10]. Queuosine is widely distributed in bacteria, and it 
is present in the first base of the anticodon of tRNA Asp , 
tRNA Asn , tRNA His and tRNA Tyr [12]; however in E. coli 
only tRNA Asp is a substrate for the GluQ-RS enzyme. 

The presence of modifications within the anticodon 
loop of the tRNA, could enhance the accuracy of 
the codon binding [13]. Then the tRNA AspQ34 might 
improve recognition of both GAC and GAU codons [14] 
and stimulate the binding of the GAU codon to the ribo- 
some [15]. In Shigella flexneri it has been shown that 
mutations in genes required for tRNA modifications, 
miaA and tgt decreased virulence. miaA is required 
for 2-methylthio-N 6 -isopentenyladenosine modification 
at position 37 of the anticodon loop and tgt is involved 
in queuosine modification at position 34 within the anti- 
codon loop [16-18]. In this study, we determined the 
role of the genome organization and its effect on the 
expression of the gluQ-rs gene in the major human 
pathogen, S. flexneri. 

Results 

Genomic organization of the S. flexneri gluQ-rs gene 

GluQ-RS is required for the synthesis of the modified 
nucleoside, GluQ, present on tRNA Asp [10,11]. By 
searching the bacterial protein database Uniprot (http:// 
www.uniprot.org/), we were able to identify GluQ-RS in 
more than a hundred bacterial species, primarily proteo- 
bacteria (Figure 1, filled symbols). From the phylogenetic 
analysis we can distinguished the three subgroups of 
enzymes described by Dubois et al, 2004 [11], which are 
characterized by the presence of the signature HXGS, 
HXGN or HXGH in the adenylate binding site. A similar 
tree was obtained using the Neighbor joining method. 



Phylogenetic analysis within the subgroup of enzymes 
with the HXGN motif, included representatives from 
the Firmicutes bacterial group (Figure 1, open square) 
together with Desulfovibrio vulgaris and Truepera radio- 
victrix enzymes. From the alignment, these members have 
8 characteristic amino acids, G 70 PDXGGXX, that do not 
align with the other GluQ-RS (Figure 1, numbering is 
derived from D. vulgaris enzyme). Further genomic ana- 
lysis indicated that the gluQ-rs gene is found primarily in 
two genomic arrangements, either alone or located imme- 
diately downstream of dksA. Searching within the String 
database [19] and GenomeNet [20], we found that the 
dksA-gluQ-rs gene organization was conserved in more 
than 40 different species, all of which were within the 
gammaproteobacteria group. These included species of 
Aerornonadales, Alter omonadales, Enterobacteriaceae, in- 
cluding E. coli and S. flexneri, Pseudomanadales, and 
Vibrionaceae (Figure 1). 

A bioinformatics analysis of the intergenic region 
between dksA and gluQ-rs showed great variation in the 
distance between the two genes among these bacterial 
species. In S. flexneri the intergenic region between the 
stop codon of dksA and the first codon of gluQ-rs is only 
39 base pairs. Therefore, we suspected that the tran- 
scription of gluQ-rs was regulated by the previously 
characterized dksA promoter [22]. To test this hypoth- 
esis, we isolated total mRNA and performed RT-PCR to 
identify an mRNA that included both genes (Figure 2A). 
The observation that there is an mRNA species contain- 
ing both genes (Figure 2 A, lane 1) indicates that they are 
co-transcribed and that the expression of gluQ-rs may 
be regulated by the dksA promoter. 

S. flexneri gluQ-rs gene is co-transcribed with dksA gene 

Although S. flexneri gluQ-rs can be transcribed from the 
dksA promoter, this did not rule out the presence of an 
additional, independent promoter. Therefore, the expres- 
sion of each gene was measured by RT-PCR during dif- 
ferent stages of S. flexneri growth in Luria Bertani (LB) 
at pH 7.4. The analysis of the dksA and gluQ-rs tran- 
scripts shows that for both mRNAs, the level is stable 
during the growth curve, with an increase of 1.3-fold at 
stationary phase compared to the early mid log phase 
(Figure 2B, compare lanes 2 and 4). However, the mRNA 
that includes the intergenic region showed variation de- 
pending on the stage of growth, increasing 20-fold at 
stationary phase compared with its expression at early 
mid log phase (Figure 2B dksA/gluQ-rs, compare lanes 2 
and 4). In order to confirm those results, a transcrip- 
tional fusion strategy was used. Different segments of 
the operon were cloned and fused to the lacZ reporter 
gene in pQF50, and promoter activity was assayed by 
(3-galactosidase activity [23]. Kang and Craig, 1990 [22] 
identified three promoters for dksA. By mean of 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 GluQ-RS is distributed within the bacterial domain. Rooted Phylogenetic analysis of selected sequences of GluQ-RS, showing the 
presence of this enzyme in the bacterial domain. Searching within the Uniprot database (http://www.uniprot.org/) of homologues protein to 
GluQ-RS were identified and then were searched in the GenomeNet (http://www.genome.jp/) database for confirmation and analysis of the 
genomic organization. Bootstrap values (>50%) where calculated by 400 replicates using the maximum-likelihood methods in the MEGA5 
software [21] and rooted with archaeal GluRS from Methonocoldococcus jonnoschii and Archoeoglobus fulgidus (not shown). In yellow background 
are shown bacterial species (in red and underlined) that are representatives having the genomic organization of dksA-gluQ-rs genes. The 
signature of each subgroup identified previously [11] is indicated. Filled symbols representing proteobacteria groups, open symbols represent 
other bacterial groups. ■: alphaproteobacteria, ▲: betaproteobacteria, •: gammaproteobacteria, ♦: deltaproteobacteria, O: actinobacteria, 
A: cyanobacteria, □: firmicutes, 0: others. 

V ) 



bioinformatics tools, including BPROM from the Soft- 
berry software package (http://linuxl.softberry.com/ 
berry.phtml), we identified those promoters in S. flexneri 
and included all three promoters in the constructs indi- 
cated in Figure 3A. The plasmid containing a fragment 
from the dksA promoters to the beginning of the gluQ-rs 
gene, with the first five amino acids of GluQ-RS, named 
pVCPDT, represents the full length dksA gene with its 
native promoters (Table 1 and Figure 3A). A second 
fusion construct, pVCDT, contains sequence from the 
beginning of the coding region of dksA through the 
beginning of gluQ-rs and also included the first five 
amino acids of GluQ-RS. Because pVCDT does not have 
the dksA promoter region, it served as the reporter for 
transcription from an independent gluQ-rs promoter. A 
third construct, pVCPD, contained the segment from 
the dksA promoter to the end of the dksA gene, hence 
this plasmid does not have the intergenic region, nor the 
first amino acids of GluQ-RS (Table 1). Each of the 
recombinant plasmids was transformed into S. flexneri 
and the p-galactosidase activity was measured when the 
bacterial cells reached mid-log phase. Analysis of the 
enzymatic activity of these reporter fusions showed that 
the strain carrying pVCDT had baseline levels of the 
enzyme (Figure 3B), indicating that there is not an inde- 
pendent promoter for gluQ-rs. Thus, the promoter 
upstream of dksA is responsible for the expression of 
both genes, at least under the conditions assayed (see 
lane pVCPDT Figure 3B). Therefore, these two results 
(Figure 2 and Figure 3B) indicate that dksA and gluQ-rs 
form an operon, and gluQ-rs lacks an additional, separ- 
ate promoter. A surprising observation was obtained 
with the clone containing pVCPD, which showed a ten- 
fold increase in enzymatic activity compared to pVCPDT 
(Figure 3B). This suggested the presence of a terminator 
or other regulatory sequence in the intergenic region 
that modulated the expression of gluQ-rs. 

The S. flexneri gluQ-rs gene has an upstream transcription 
terminator 

In order to explain the difference observed in expression 
of lacZ from the recombinant plasmids pVCPDT and 



pVCPD a bioinformatic analysis using mFold [26] was 
performed to search for possible secondary structures in 
the mRNA. A potential transcriptional terminator was 
found at the beginning of the gluQ-rs gene, leaving the 
first predicted AUG codon located on the bulge of this 
terminator (Figure 4A). In order to determine the func- 
tionality of this terminator, we performed site directed 
mutagenesis to disrupt the structure in the predicted 
stem (Figure 4A). As shown in Figure 4B, the plasmid 
containing the mutations, pVCPDTMut had >2-fold 
higher enzymatic activity (^><0.05) than the plasmid 
containing the wild type sequence. This result suggested 
that the intergenic region upstream of gluQ-rs contains 
a transcriptional terminator. 



Identification of the first methionine 

The first methionine in the predicted GluQ-RS protein 
corresponds to the one located on the bulge of the ter- 
minator structure (Figure 4A), which also contains a 
possible Shine-Dalgarno sequence. However, in related 
species like Escherichia fergusonii that also have the 
terminator structure, a methionine is not present at that 
location. In the S. flexneri sequence, there is another 
AUG codon in the same reading frame 27 nucleotides 
downstream from the one in the terminator. In order to 
determine which methionine is the start site for transla- 
tion of the S. flexneri GluQ-RS, we constructed a vector 
that included the intergenic region from the stop codon 
of the dksA gene to the end of gluQ-rs cloned into the 
expression vector pET15c. This allowed expression of 
C-terminal His-tagged GluQ-RS under T7 promoter 
control. The protein was partially purified by affinity 
chromatography as described elsewhere [10], and the 
sequence of the amino terminus of the GluQ-RS 
protein was determined to be NH 2 -T-D-T-Q-Y-I-G-R- 
F-A-P, which corresponds to the amino acid sequence 
after the second methionine. Therefore, the initiator 
methionine is not the one indicated in the database, 
and the protein is 298 amino acids. Surprisingly, there 
is no obvious Shine-Dalgarno sequence adjacent to the 
initiator methionine we identified (Figure 5). 
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Figure 2 gluQ-rs is co-transcribed with dksA in S. flexneri 2457T. A) Agarose gel showing the amplified product of the full-length operon 
extending from the dksA gene through the end of g I 'uQ-rs (1.44 kpb). Total RNA isolated during mid log phase growth of S. flexneri was subjected 
to reverse transcriptase and PCR (RT-PCR) using primers opeF/opeR (Table 2). M: molecular marker, sizes are indicated. Lane 1: RNA treated with 
reverse transcriptase, Lane 2: genomic DNA isolated from S. flexneri 2457T, Lane 3: RNA without reverse transcriptase treatment, Lane 4: negative 
control of PCR reaction without DNA. B) Electrophoretic analysis of each amplified gene fragment, dksA (dksAF/dksAR; 436 bp), gluQ-rs (gQRSF/ 
gQRSR; 508 bp), the intergenic region dksA/gluQ-rs (interF/interR; 496 bp) and the ribosomal RNA 16S {rrsHF/rrsHR, 589 bp). Total RNA isolated 
during different phases of the growth curve of S. flexneri 2457T was subjected to RT-PCR to detect the corresponding fragment. Lane V. lag 
phase, Lane 2: early mid log phase, Lane 3: mid log phase, Lane 4: stationary phase. +: corresponds to amplification using genomic DNA. RNA: 
Isolated RNA without reverse transcriptase treatment. -: negative control PCR reaction without DNA. Each band was estimated using Image J 
software (V1.46) and its amount was estimated using the fragment corresponding to 500 bp of the DNA ladder. 



Phenotype of the S. flexneri gluQ-rs mutant 

To determine the role of GluQ-RS in S. flexneri growth 
and virulence, a deletion mutant of the gluQ-rs gene was 
constructed in S. flexneri 2457T. The mutant was com- 
pared to the wild type by Biolog phenotype MicroArrays 
(Biolog, Inc., Almeda, CA). The major difference 
observed for the mutant was impaired metabolism when 
grown under osmotic stress conditions (Figure 6). The 
mutant had a longer lag and reduced growth compared 
to the wild type in the presence of increasing concentra- 
tions of potassium chloride, sodium sulfate, sodium 
formate, sodium benzoate, sodium nitrate and sodium 
nitrite. The phenotype was complemented with the gluQ- 
rs gene cloned into an expression vector. No differences 
were observed in the growth or metabolism of these 
strains when they were incubated in presence of 1% 
sodium chloride, which was similar to LB (Figure 6 and 
data not shown). 

Because expression of dksA is required for S. flexneri 
virulence [27], and growth of Shigella in the intracellular 
environment may induce a stress response, we also mea- 
sured invasion and plaque formation by the gluQ-rs 
mutant. However, no significant differences were noted 
(data not shown), suggesting that GluQ-RS is not essen- 
tial for invasion or intracellular growth of S. flexneri. 



Discussion 

Conserved dksA-gluQ-rs genomic organization 
in gammaproteobacteria 

GluQ-RS, a paralog of GluRS synthetase, is involved in 
the formation of GluQ, the nucleoside located at the wob- 
ble position of tRNA Asp in bacteria. The protein is present 
in Firmicutes, Actinobacteria, Cyanobacteria, Alphapro- 
teobacteria, Betaproteobacteria, Gammaproteobacteria 
and Deltaproteobacteria (Figure 1). From the phylogenetic 
analysis we distinguished the three subgroups described 
previously [11] based on the HIGH motif that is present 
in the class I aminoacyl-tRNA synthetases [2]. As was 
described previously [11], all GluQ-RS enzymes are charac- 
terized by the replacement of a threonine in GluRS 
enzymes, which is involved in the recognition of the amino 
acid and the terminal adenosine of the tRNA Glu (Thr 133 of 
Methanocaldococcus jannaschii GluRS enzyme) by isoleu- 
cine, leucine or valine at that position (Ile 47 of S. flexneri 
GluQ-RS). This substitution is also conserved in all 
enzymes analyzed here, including those from the Firmicutes 
group. 

The gluQ-rs gene is widely distributed in the bacterial do- 
main; however, its genome organization is variable. We 
observed that only in members of the gammaproteobacteria, 
namely Aerornonadales, Alteromonadales, Pseucfamonadaceae, 
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Figure 3 The transcription of the gluQ-rs gene is controlled by a termination stem loop. A) Schematic representation of the operon, the 
arrows indicate the position of each promoter identified by our bioinformatics analysis and experimentally determined by Kang and Craig, 1990 
[22]. The putative p-independent terminator is represented by the stem loop symbol upstream of gluQ-rs gene. The horizontal bar represents 
the DNA region amplified and cloned into pQF50 (Table 1). The recombinant plasmids are described in Table 1. pVCDT does not have the dksA 
promoter but has the terminator. pVCPDT has the promoter region of dksA and the terminator upstream of g 1 'uQ-rs; therefore, it represents the 
genomic organization of the operon. pVCPD also has the promoter of dksA but lacks the terminator region. The size of each fragment is 
indicated. B) (3-galactosidase activity of each protein extract obtained from the corresponding clone. The data represent the average of three 
experiments in triplicates and the Student t test was used to compare the means between each clone with the empty vector. *** p values <0.05 
were considered statistically significant. 



Enterobacteriaceae and Vibrionaceae, the gluQ-rs gene is 
located immediately downstream of the dksA gene (Figure 1). 
A more detailed analysis shows that even within this genomic 
organization there are differences. In some species of Pseudo- 
monadaceae, such as P. aeruginosa, P. entomophila, and P. 
fluorescens, we observed the same genomic structure as in E 
coli or S. jlexneri, with a distinctive terminator between the 
genes. In contrast, while the dksA gene is also upstream of 
gluQ-rs in some P. syringae, there are insertions of an encoded 
transposase or more than a 400 base pairs separating both 
genes without a detectable terminator. However, using bio- 
informatics tools we detected a possible promoter within this 
region in P. syringae (data not shown), indicating that the 
expression of the gluQ-rs gene may be under control of its 
own promoter, a question that remains to be addressed. 

Stringent response and tRNA modification 

Our bioinformatic analysis shows the presence of a tran- 
scriptional terminator and lack of a g/wQ-rs-specific 



promoter. This is consistent with our results, in which we 
did not detect any activity from promoters other than 
those upstream of the dksA gene (Figure 3). This unusual 
arrangement suggests that gluQ-rs expression is dependent 
on d/<sA-regulated conditions. Because DksA is a key 
member of the stringent response in bacteria and regulates 
a number of processes in the cell, including its own 
expression [25,28], the data suggest that there is coordinate 
regulation of tRNA modification and other DksA targets. 
Although we could not detect any promoter activity spe- 
cific for gluQ-rs in the growth conditions tested (i.e. alter- 
ing the pH, presence of glutamate), we cannot discount 
the possibility that the gene is specifically regulated under 
some other conditions. The regulon database (http://regu- 
londb.ccg.unam.mx/index.jsp) indicates that the E. coli 
gluQ-rs gene has a recognition site for the a 24 subunit of 
RNA polymerase. From our analysis, this sequence is iden- 
tical to S. flexneri, but there is no experimental evidence of 
this recognition. Interestingly, when the gluQ-rs gene was 
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Table 1 Bacterial strains and plasmids used in this work 

Bacterial strains or plasmid Characteristics 

Shigella flexneri 

S. flexneri 2457T Wild type strain 

S. flexneri 2457T AgluQ-rs::kan Deletion mutant of gluQ-rs gene 

Escherichia coll 

E. coll W31 1 0 AgluQ-rs::kan Deletion mutant of gluQ-rs gene 



deleted in S. flexneri, the mutant showed impaired growth 
in the presence of osmolytes (Figure 6). A recent publica- 
tion demonstrated that a 24 and cr s proteins from Salmon- 
ella enterica serovar Typhi are important for the 
expression of several genes induced by osmotic stress in 
this bacterium [29]. Moreover, the expression of the gene 
encoding a 24 in E. coli is regulated by the stringent 
response [30]. The possible role of o* 24 on the expression 
of gluQ-rs under osmotic stress might be interesting 
to study. 

GluQ-RS is an enzyme responsible for the formation 
of the GluQ tRNA modification, and two independent 
groups [10,11] have shown that this enzyme required a 
high concentration of glutamate to be activated and 
transferred to the queuosine base present on the 
tRNA Asp . Interestingly, one of the first events to occur 
when bacteria are subject to high osmolyte stress is an 
increase in glutamate levels within the cytoplasm [31]. 
Our observation indicates an important role of the tRNA 
modification for the growth of S. flexneri in the presence 
of osmolytes (Figure 6). Other tRNA modifications 
might play a similar role in this stress condition. In E. 



Source or reference 



Laboratory stock 
This work 

[10] 



coli, inactivation of the yfiC gene, responsible for the 
modification at the adenosine 37 present on the tRNA Val , 
leads to a high sensitivity to osmotic stress [32]. 

Transcription of gluQ-rs is regulated by a terminator 

The results obtained in the present work show the pres- 
ence of a terminator and suggested the functionality of 
this structure (Figure 3 and Figure 4). To our knowledge, 
there are few examples of bacterial genes that have simi- 
lar structures. There is a terminator structure upstream 
of the DNA primase gene, dnaG, which also has an un- 
usual Shine Dalgarno sequence [33]. Another example is 
the recX gene in E. coli, where readthrough accounts for 
approximately 10% of its transcription [34]. 

The two characteristics of gluQ-rs described in this 
work, co-transcription with the upstream gene and the 
presence of a terminator immediately upstream, allow us 
to propose that both the transcription and translation 
process could be regulated in the gluQ-rs gene. It has 
been described, that the presence of terminators up- 
stream of the coding region might be part of a regula- 
tory system such as a riboswitch [35]. Riboswitches for 
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F cp80lacZAMI5 A(lacZYA-argF) UI69 recA I endA I hsdR 1 7 (rK-, mK+) [24] 
phoA supE44 X-thi-l gyrA relAI 

F ompTgal dcm Ion hsdS B (r B m B ) X(DE3 [lad lacUV5-T7 gene I indl Invitrogen 
sam7 nin5]) 



bla, pMB1 ori, lacZ peptide, fl phage ori Fermentas 

bla, pMB1 ori, lacZ gene without promoter [23] 

Empty vector, a modified version of pET15b This work 

S. flexneri fragment from nucleotide +58 a of dksA gene to beginning This work 

of g 1 'uQ-rs gene (+590) cloned into pQF50. Pair of primers used were 

PgluQF/PdksARCT. 

5. flexneri fragment from nucleotide -506 of dksA gene to beginning This work 

of gluQ-rs gene (+590) cloned into pQF50. Pair of primers used were 

PdksAF/PdksARCT. 

S. flexneri fragment from nucleotide -506 of dksA gene to beginning This work 

of g 1 'uQ-rs gene (+590) cloned into pQF50, with the terminator mutated 
by the nucleotides indicated in Figure 4a. 

S. flexneri fragment from nucleotide -506 of dksA gene to This work 

nucleotide +527 (end of dksA gene) cloned into pQF50. Pair of primers used 
were PdksAF/PdksARST. 

5. flexneri gene from nucleotide +509 (stop codon of dksA) to This work 

nucleotide +1469 (last codon of gluQ-rs without stop codon). Pair of primers used 
were ATGGQRSF/ATGGQRSR. 



fragments cloned are numbered based on the transcription start of dksA identified by [25]. 
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Figure 4 Functionality of the transcriptional terminator upstream of gluQ-rs. A) Schematic representation of the terminator with a 
AG = -14.7 Kcal/mol identified using Mfold software [26]. Bases shaded in grey indicate the two possible AUG start codons, one located in the 
bulge of the terminator structure and the other located 27 nucleotides downstream. The arrows indicate the site directed mutagenesis location, 
with the corresponding nucleotide changes designed to disrupt the predicted structure. B) (3-galactosidase activity of protein extracts obtained 
from the corresponding clones. The plasmid pVCPDTMut has a similar construction as pVCPDT but contains the mutated terminator indicated 
above. The data represent the average of three experiments, each done in triplicate, and the Student t test was used to compare means between 
the pVCPDT and pVCPDTMut clones. *** p values <0.05 were considered statistically significant. 



genes involved in queuosine formation have been 
described, in which the precursor preQl is the ligand of 
the mRNA structure [36]. Using the riboswitch server 
(ribosw.mbc.nctu.edu.tw [37]), we did not identify any 
potential riboswitch (data not shown). However we can- 
not discount that the terminator described here might 
be part of a regulatory circuit similar to a riboswitch, 
or that an unidentified protein might bind the termin- 
ator structure. 

GluQ modification and codon bias 

tRNA modifications present at the anticodon loop might 
be important for the accuracy of codon reading during 
the translation processes [13]. Morris et al, 1999 [14] 
proposed, based on molecular modeling, that the 
tRNA AspQ34 might improve recognition of both GAC 
and GAU codons, consequently the interaction of the 
codon GAU with the anticodon of tRNA AspG34 could be 
less efficient. In fact, in S. flexneri there are a few genes 



such as sit A j virF and proX (an inducible gene under 
osmotic stress) that have a bias toward those codons 
that favor the modified tRNA. Thus, while there is no 
obvious loss of plaque formation in the gluQ-rs mutant 
(data not shown), the absence of GluQ-RS may influence 
the expression of proteins such as SitA that are required 
for fitness of Shigella in the host [38]. 



Conclusions 

In this work we have shown that the expression of gluQ-rs, 
a gene codifying an enzyme involved in the formation of 
GluQ present on the tRNA Asp , is under the control of the 
dksA promoter. Also, we show the presence of a func- 
tional terminator that controls the expression of gluQ-rs. 
Finally, we present data that suggest that the presence of 
modification of the tRNA Asp is important for survival of 
the human pathogen Shigella flexneri under osmotic stress 
conditions. 
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Figure 5 Determination of the first methionine of GluQ-RS. The cloning strategy utilized is shown at the top. A fragment from the stop 
codon of dksA to the end of gluQ-rs gene was amplified from S. flexneri genomic DNA with the primers ATGGQRSF/ATGGQRSR and cloned 
into a pET15c vector using the restriction sites BomVW and Xho\. Therefore this clone represents the operon, but dksA was replaced with the 
plasmid encoded fragment pelB. The transcription of this plasmid, named pATGGQRS, is controlled by the T7 promoter and translation is 
controlled by the Shine-Dalgarno (SD) sequence of pel % both contained within the plasmid. The GluQ-RS protein synthesized has a histidine 
tag (H 6 ) at the C-terminus, facilitating protein purification. The putative p-independent terminator is represented by the stem loop symbol 
upstream of gluQ-rs gene , and the black box in the gene indicates the two possible peptides, depending of which methionine is utilized. 
Bottom: The SDS -polacrylamide gel electrophoresis showing the supernatant extract (S) and the partially purified protein (Pur) produced by 
cells carrying the recombinant plasmid. The predominant band indicated by the arrow was excised and subjected to amino terminal sequencing, 
yielding the following amino acid sequence: T-D-T-Q-Y-l-G-R-F-A-P. This corresponds to the sequence following the second methionine. The 
size of the molecular markers (M) are given in kDa. 



Methods 

Bacterial growth conditions 

The bacterial strains and plasmids used in this study are 
described in Table 1. E. coli strains were maintained on LB- 
agar (10 g of tryptone per liter, 5 g of yeast extract per liter, 
10 g of NaCl per liter and 15 g of agar per liter), Shigella 
strains were maintained on Trypticase Soy Agar plus 0.01% 
congo red. All strains were stored at -80°C in LB broth plus 
20% glycerol. The bacteria were grown in LB broth adjusted 



to pH 7.4 with 40 mM MOPS (Merck) or M9 minimal 
media [24], When necessary, ampicillin was added to a final 
concentration of 100 ug/ml. Bacterial growth was moni- 
tored by optical density at 600 nm (OD 600 ). 

Bioinformatics tools to construct the phylogenetic tree 

The protein sequences were obtained from the Uniprot 
database (http://www.uniprot.org/) and then were searched 
in the GenomeNet (http://www.genome.jp/) to confirm the 
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Figure 6 The gluQ-rs mutant is sensitive to growth in high osmolarity. Biolog phenotypic MicroArrays were used to characterize the growth 
and metabolism of the gluQ-rs mutant and its wild type parent, 2457T. Wild type (black line) transformed with the empty plasmid, 2457T Ag/uQ- 
rsr.kan (red line) transformed with the empty plasmid pCM and 5. flexneri AgluQ-rs::kon (blue line) transformed with the complementary plasmid 
pTRCGQ were incubated for 24 h with medium in the corresponding plate for osmolytes (PM9). A) sodium chloride 1% B) sodium benzoate 
20 mM pH 5.2. C) sodium nitrate 100 mM. Metabolism was monitored by measuring reduction of the tetrazolium dye in the medium at 15 min 
intervals and is shown as units. 
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genomic organization. A selected number of GluQ-RS 
enzymes were aligned using the MUSCLE algorithm [39] 
and analyzed using the maximum-likelihood method 
based on the JTT matrix-based model. The percentage 
of trees in which the associated proteins clustered to- 
gether is shown next to the branches. The analysis 
involved 54 amino acid sequences, including the GluRS 
proteins from Methanocaldococcus jannaschii and 
Archaeoglobus fulgidus as an outgroup. All positions 
containing gaps and missing data were eliminated. There 
were a total of 199 positions in the final dataset. Evolu- 
tionary analyses were conducted in MEGA5 [21]. 



RNA isolation and synthesis of cDNA 

Total mRNA was obtained during the growth of S. flexneri 
2457T using the RNeasy mini kit following the supplier 
instructions (Qiagen). The purified nucleic acid was trea- 
ted with RNase- free DNase (Fermentas) and its concen- 
tration was estimated by measuring the optical density at 
260 nm (OD 2 6o)« Approximately 1 ug of total RNA was 
subjected to reverse transcription using M-MuLV poly- 
merase (Fermentas) and random primers following the 



providers protocol. The cDNA was amplified using 
specific PCR primers for each gene of interest (Table 2). 

Construction of transcriptional fusions 

Transcriptional fusions were constructed to study the ex- 
pression control of gluQ-rs. Fragments of the dksA-gluQ-rs 
region were fused to lacZ in the vector pQF50 by using the 
BamHl and HindLll restriction sites [23]. Each fragment 
was amplified from S. flexneri genomic DNA using the 
indicated primers (Tables 1 and 2) with the High Fidelity 
PCR Enzyme Mix polymerase (Fermentas) and cloned into 
pQF50 (Table 1). Once the sequence of each clone was 
confirmed, the recombinant plasmid was introduced into S. 
flexneri 2457T by electroporation. The nomenclature of the 
recombinants plasmids is: P for promoter of the dksA gene, 
D for the dksA gene and T for a terminator structure. 

P-galactosidase activity 

S. flexneri transformed with the corresponding con- 
structs were cultured overnight in LB, a 1:50 dilution 
was inoculated into 10 ml culture of LB pH 7.4 and 
grown to an OD 600 of 0.5. Aliquots of 0.5 ml of each 
strain containing the clone or the empty vector were 



Table 2 Primer sequences 



Name 



Sequence 5'- 3 /a 



Reference and characteristics 



opeF TAAGGAGAAGCAACATGCAAGA 

opeR ATAGCTCAGCATGACGCAm 

dksAF ATGCAAGAAGGGCAAAACCG 

dksAR GCGAAmCAGCCAGCGm 

interF AGTGGAAGACGAAGAmCG 

i nterR TCOTGTOATGTAACC AGG 

gQRSF TTCAAAGAGATGACAGACACACAG 

gQRSR CACGGCGATGAATGATAAAATC 

rrsHF CCTACGGGAGGCAGCAG 

rrsHR CCCCCGTCAATOCmG AGm 

pen BR GATGGAGCCGAAAATGTOT 

PdksAF GGATCCAAGCGAAGTAAAATACGG 

Pd ksARST AAGCTTGTG ATGGAACGGCTGTAAT 

Pd ksARCT AAGCTTCTGTGTGTCTGTCATCTCmG 

PgluQF GGATCCAAGAAGGGCAAAACCGTA 

TERGQ2 CCTTAI I I I I I GTOAAAGAGATGACAGACACACAGA 

TERMGQ3 ATAAGGCGGGAGCATAAC GGAGGA GTGGTAAAC 

M1 3R GCGGATAACAAmCACACAGG 

ATGGQRSF GGATCCGTAATOCAGCCGTCCATC 

ATGGQRSR CTCGAGGCATGACGCA^TGAGAATG 

virFF AGCTCAGGCAATGAAACmGAC 

virFR TGGGOTGATATOCGATAAGTC 



This work. RT-PCR of dks A operon from nucleotide +40 to +1477 b 

This work. RT-PCR of dksA gene from nucleotide +54 to +488 

This work, RT-PCR of intergenic region from nucleotide +368 to +863 

This work, RT-PCR of g I uQ-rs gene from nucleotide +567 to +1074 

[40] RT-PCR of ribosomal RNA 16S 

Reverse of pen B gene from nucleotide +1993 
BamVW site, from nucleotide -506 
HindlW site, to nucleotide +527 
HM\\\ site, to nucleotide +590 
BamVW site, from nucleotide +58 
Recognition from nucleotide +555 

Recognition from nucleotide +560, underline sequence are nucleotides changed 
Recognition site in pTZ57R^~ 

BomVW site, from nucleotide +507. Underline nucleotides correspond to the stop 
cod on of dks A 

Xho\ site, to nucleotide +1469 
[41] 



Nucleotides in bold are indicated restriction site. 

fragments cloned are indicated based on the transcription start of dksA identified by [25]. 



Caballero et al. BMC Microbiology 2012, 12:226 
http://www.biomedcentral.com/1471 -21 80/1 2/226 



Page 11 of 13 



assayed for p-galactosidase activity according to Miller 
[42]. The data were analyzed using the software Graph- 
Pad Prism V5.01. 

Site directed mutagenesis 

A possible transcription terminator between dksA and 
gluQ-rs was identified using the program Mfold [26]. Site 
directed mutagenesis by overlap PCR was performed to 
disrupt the predicted terminator [43]. Using the fragment 
VCPDT cloned in the vector pTZ57R/T as template, was 
amplified a 1,072 bp fragment, which include the mutation, 
using the primers PdksAF and TERMGQ3, while a second 
fragment of 162 bp overlapping the mutated region, was 
obtained with primers TERGQ2 and M13R (Table 2). Both 
fragments (1,072 bp and 162 bp) were digested with Dpnl, 
purified and mixed at equimolar quantities to carry out a 
PCR reaction using the 5' and 3' ends primers (PdksAF and 
PdksARCT). The 1,110 bp amplified fragment was 
cloned in the vector pTZ57R/T and sequenced to 
verify the mutation. This plasmid was digested with 
BamHl and Hindlll and the fragment subcloned in to 
the vector pQF50. 

Determination of first methionine of GluQ-RS 

In order to establish which is the first AUG codon of 
gluQ-rs, the recombinant plasmid pATGGQRS was 
constructed. A PCR reaction was performed using the 
primers ATGGQRSF and ATGGQRSR (Table 2) and 
genomic DNA from S. flexneri. The amplified fragment, 
containing the BamHl site, stop codon of dksA, the 
intergenic region with the terminator, the gluQ-rs read- 
ing frame without its stop codon and the Xhol site was 
cloned into pET15c, a modified version of pET15b, 
which was constructed by inserting the 290 bp Xbal and 
Blpl fragment of pET20b containing the polylinker into 
pET15b. This construct allowed the synthesis of a C- 
terminal histidine tagged protein under the transcription 
control of the T7 promoter. The construct was trans- 
formed in BL21(DE3) strain and the His-tagged protein 
was partially purified by affinity chromatography as 
described previously [10]. The eluted protein was trans- 
ferred to a PVDF membrane and stained with Coomassie 
blue. The predominant band of the expected size 
(34.6 kDa) was sequenced at the Protein Core Facility of 
the Institute for Cellular and Molecular Biology, University 
of Texas at Austin. 

Construction of the plasmid for complementation 
of the gluQ-rs mutation 

This plasmid was constructed from the pATGGQRS 
plasmid in which the T7 promoter was removed by 
digestion with Bglll and Ncol enzymes and replaced by 
the TRC promoter obtained from pTRC99a plasmid 
by amplification and digestion with BamHl and Ncol to 



obtain the pTRCGQ plasmid. The empty plasmid (pCM) 
was constructed by incorporating the TRC promoter 
into the pET15c plasmid. 

Inactivation of gluQ-rs gene in S. flexneri 

Deletion of gluQ-rs was carried out using the A red 
recombinase method [44] with the following modifica- 
tions. S. flexneri 2457T carrying pKD46 and prepared as 
described elsewhere [44] was transformed with a purified 
PCR fragment amplified from the E. coli AgluQ-rs::kan 
mutant strain using primers dksAF and pcnBR (Table 2), 
increasing the homologous DNA region to more than 
450 bp at each side. The mutant was isolated following 
overnight growth at 37°C on LB-agar containing kana- 
mycin (50 (ig/ml). The deletion was confirmed by PCR 
using the same pair of primers (dksAF-pcnBR) and using 
each primer together with an internal primer as 
described previously [44] . The presence of the S. flexneri 
virulence plasmid was also confirmed by PCR amplifica- 
tion of the virF gene using primers virFF and virFR 
(Table 2). 

Effect of the absence of gluQ-rs gene in S. flexneri 
metabolism 

The effect of the deletion of the gluQ-rs gene on the me- 
tabolism of S. flexneri was analyzed by Biolog phenotype 
MicroArrays following the manufacturers instructions 
(Biolog, Inc., Almeda, CA). Strains were grown at 30°C 
overnight and 5 ml of LB was inoculated with a 1:100 dilu- 
tion and grown at 37°C to reach an OD 650nm of 0.5. The 
cells were then washed and resuspended to 2.5 x 10 7 cfu/ 
ml and diluted 200 fold in to a solution of IF- 10a medium 
(Biolog). Each well was inoculated with 1.2 x 10 4 cfu 
(0.1 ml per well) into the corresponding plates and incu- 
bated for 24 hrs at 37°C. The metabolism was recorded 
and analyzed by the Omnilog software (V 1.20.02) (Biolog, 
Inc., Almeda, CA). 
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